Quantum-optical communication rates through an amplifying random medium 
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We study the competing effects of stimulated and spontaneous emission on the information 
capacity of an amplifying disordered waveguide. At the laser threshold the capacity reaches a 
"universal" limit, independent of the degree of disorder. Whether or not this limit is larger or 
smaller than the capacity without amplification depends on the disorder, as well as on the input 
power. Explicit expressions are obtained for heterodyne detection of coherent states, and generalized 
for arbitrary detection scheme. 
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To faithfully transmit information through a commu- 
nication channel, the rate of transmission should be less 
than the capacity of the channel 0,0] . Although current 
technology is still far from the quantum limit, there is an 
active scientific interest in the fundamental limitations to 
the capacity imposed by quantum mechanics ||||]. Ulti- 
mately, these limitations originate from the uncertainty 
principle, which is the source of noise that remains when 
all external sources have been eliminated. 

An important line of investigation deals with strategies 
to increase the capacity. One remarkable finding of recent 
years has been the beneficial role of multiple scattering by 
disorder, which under some circumstances can increase 
the capacity by increasing the number of modes (scatter- 
ing channels) that effectively carry the information . 
Quite generally, the capacity increases with increasing 
signal-to-noise ratio, so that amplification of the signal is 
a practical way to increase the capacity. When consider- 
ing the quantum limits, however, one should include not 
only the amplification of the signal (e.g. by stimulated 
emission), but also the excess noise (e.g. due to sponta- 
neous emission). The two are linked at a fundamental 
level by the fluctuation-dissipation theorem, which con- 
strains the beneficial effect of amplification on the capac- 

ity §■ 

While the effects of disorder and amplification on com- 
munication rates have been considered separately in the 
past, their combined effects are still an open problem. 
Even the basic question "does the capacity go up or down 
with increasing gain?" has not been answered. We were 
motivated to look into this problem by the recent inter- 
est in so called "random lasers" ||||. These are optical 
media with gain, in which the feedback is provided by 
disorder instead of by mirrors (as it is in a conventional 
laser). Below the laser threshold these materials behave 
like linear amplifiers with strong intermode scattering, 
and this results in some unusual noise properties fToyilf , 
As we will show here, the techniques developed in con- 
nection with random lasers can be used to predict under 
what circumstances the capacity is increased by amplifi- 
cation. 

We consider the transmission of information through a 



linear amplifier consisting of an iV-mode waveguide that 
is pumped uniformly over a length L (see Fig. [I]). We will 
refer to amplification by stimulated emission, but one can 
equally well assume other gain mechanisms (for example, 
stimulated Raman scattering The amplification oc- 

curs at a rate l/r a . The waveguide also contains passive 
scatterers, with a transport mean free path I. The com- 
bined effects of scattering and amplification are described 
by a 2N x 2N scattering matrix S which is super-unitary 
(SS' — 1 positive definite). 

We assume that the information itself is of a classical 
nature (without entanglement of subsequent inputs), but 
fully account for the quantum nature of the electromag- 
netic field that carries the information. The quantized 
radiation is described by a vector a ln of bosonic annihila- 
tion operators for the incoming modes and a vector a out 
for the outgoing modes. The two vectors are related by 
the input-output relation @p|Jl3] 
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The vector of bosonic creation operators describes 
spontaneous emission by the amplifying medium. The 
fluctuation-dissipation theorem relates the matrix U to 
Shy 



UU* = ss^ - a. 



(2) 



The first communication channel that we examine is 
heterodyne detection of coherent states ||. The sender 
uses a single narrow-band mode a (with frequency ojq 
and band width Aw), to transmit a complex number \x by 
means of a coherent state (such that a™|/i) = mIa*))- 
The receiver measures a complex number v by means of 
heterodyne detection of mode (3. Two sources of noise 
may cause v to differ from \x: Non-orthogonality of the 
two coherent states and \v)\ and spontaneous emis- 
sion by the amplifying medium. 

The a priori probability p(n) that the sender transmits 
the number fj,, and the conditional probability V{v\[i) 
that the receiver detects v if fx is transmitted, determine 
the mutual information B 
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Wc have defined p{v) = J d 2 fiV(v\fi)p(fi). The chan- 
nel capacity C (measured in bits per use of the chan- 
nel) is obtained by maximizing / over the a priori distri- 
bution p(fx), under the constraint of fixed input power 
P = P J d 2 /x \uPp(fJ,) (with P = hujoAuj/2-K). As 
argued in Ref. Jig ], any randomness in the scattering 
medium that is known to the receiver but not to the 
sender can be incorporated by averaging / before maxi- 
mizing, hence 



C = max (J) 
pO) 



(4) 



The brackets (• • •) indicate an average over different po- 
sitions of the scatterers. 
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FIG. 1. Communication channel consisting of an TV- mode 
waveguide that is amplifying over a length L. Both sender 
and receiver use a single narrow-band mode (indicated by a 
plane wave). 

The calculation of the capacity is greatly simplified by 
the fact that the spontaneous emission noise is a Gaus- 
sian superposition of coherent states. This is expressed 
by the density matrix of the amplifying medium, 



PmediumCX / d 2 /3 6Xp(- \0\ 2 //) \/3) (J3\ , 



(5) 



where (3 is a vector of 2N complex numbers and \0) is the 
corresponding coherent state (such that b n \(3) — (3 n \(3)). 



The variance / = N uppeT (N u 



ippcr 



lower ) 



_1 depends 

on the degree of population inversion of the upper and 
lower atomic levels that generate the stimulated emission. 
Minimal noise requires a complete population inversion: 
Mower = =>• / = 1. We will consider that case in what 
follows. 

We will similarly assume that the heterodyne detection 
adds the minimal amount of noise to the signal. (This 
requires that the image band is in the vacuum state 0.) 
The conditional probability is then given simply by a 
projection, 



V(v\n) = (v\p out (n)\v), 



of the density matrix p ont {fJ>) of the outgoing mode (3 onto 
the coherent state \v) (for an incoming coherent state \p) 
in mode a). In view of Eqs. ([!]) and @ we have (for 
/ = 1) 



PoutO) 
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W')W\- (7) 



This is again a Gaussian superposition of coherent states, 
but now the variance is related by Eq. (||) to the scatter- 
ing matrix of the medium: J2rA u f3n\ 2 = S„ \ s Pn\ 2 - 1- 
Substituting p out into Eq. (@), and using |(i/|z/)| 2 = 



exp(— \v — v'\ ), we arrive at 



V{y\p) ex exp 



Sf3a^\ 



J2n \S/3r. 



(8) 



This expression for the conditional probability has the 
same Gaussian form as in previous studies pM of com- 
munication channels degraded by Gaussian noise, but the 
essential difference is that in our case the noise strength 
is not independent of the transmitted power, but related 
to it by the fluctuation-dissipation theorem (^). 

The calculation of the capacity proceeds as in Refs. 
Jl5| , ^6| . The optimum a priori distribution p(p) oc 
exp(^|/i| 2 Po/-P) is independent of the scattering matrix 
S, so the maximization and disorder average in Eq. (^) 
may be interchanged. The result is 



c = (io g2 (i+K)), n 



{P/P*)\S Pa \ 2 



(9) 



The quantity 1Z is the signal-to-noise ratio at the re- 
ceiver's end. We can write 1Z equivalently in terms of 
the transmission matrix t (from sender to receiver) and 
the reflection matrix r (from receiver to receiver): 



K = 



(P/Po)\t 0a \ 



ELi(W 2 + M 2 ) 



(10) 



In the absence of intermode scattering one has \t Tl 



and 



= 0, hence 1Z — 8 a pP/Po and C 



log 2 (l + S a pP/Po): independent of the amount of ampli- 
fication. The increase in capacity by stimulated emission 
is canceled by the extra noise from spontaneous emis- 
sion 0| . In the absence of amplification, but in the 
presence of scattering, one has J^n 1^/3" 1 2 = 1; hence 
C = (log 2 (l + \tf3 a \ 2 P/Po)}. The capacity is reduced by 
intcr-modc scattering in the same way as for the lossy 
channel studied in Ref. |l7) . 

The average over the scatterers can be done analyt- 
ically in the limit N ^> 1 of a large number of modes 
in the waveguide. Sample-to-sample fluctuations in the 
denominator a = J2n(\t/3n\ 2 + l r ^n| 2 ) are smaller than 
the average by an order N, so these fluctuations may be 
neglected and we can replace the denominator by its av- 
erage tr. The fluctuations in the numerator r = |t/3 a | 2 
can not be ignored. These are described (for N 3> 1) by 



(g) the Rayleigh distribution V{t) = r x e T I T . Integrating 



log 2 (l + (P/Pq)t/o-) over r with distribution V(t) we 
arrive at 

C = e K <r«r(0;^- ff 1 )/ln2, TZ cS = (11) 

Poa 

with r(0; x) the incomplete gamma function. The de- 
pendence of the capacity C on the effective signal-to- 
noise ratio Tics is plotted in Fig. @. It lies always below 



2 



the capacity Cq = log 2 (l + H c s) which one would ob- 
tain by ignoring fluctuations in r. For TZ c s <C 1 the two 
capacities approach each other, C « Cq w 72- e flf/ha2, 
while for 7?. c ff ^> 1 one has Co ~ log 2 lZ c g versus 
C ~ log 2 TZcs — 7/ In 2 (with 7 w 0.58 Euler's constant). 




R eff 

FIG. 2. Capacity C for heterodyne detection of coherent 
states as a function of the signal-to- noise ratio TZ c a. The re- 
sult ([n]) lies below the value Co = log 2 (l + lZ c s) that ignores 
statistical fluctuations. 

The quantity TZ c g depends on three length scales, the 
length L of the amplifying region, the mean free path 
I, and the amplification length l a = \JDr a (with D the 
diffusion constant). The two averages f , a can be calcu- 
lated from the diffusion equation in the regime I <g. l a ,L. 
There is a weak dependence on the mode indices a, /3 in 
this diffusive regime, which we ignore for simplicity. The 
result is 

Al/3l a l-cos(Z,/i n ) 

T = AT ■ it n s i cr = 1 + (41 3l a ) , . (12) 

Nsm(L/l a ) sm(L/l a ) 

The effective signal-to-noise ratio becomes 



p 

iVPo" 



[l-cos(L/Z a ) + (3/ a /40sin(L/; a )] \ (13) 



Without amplification, for l a ^S> L, one has TZ c g — 
^(1/NL)P/Pq. Amplification increases TZ e s, up to the 
limit Ti. e ff — ► P/2NPq that is reached upon approaching 
the laser threshold l a — > L/71-. 

We conclude that amplification in a disordered wave- 
guide increases the capacity for heterodyne detection of 
coherent states, up to the limit 



Coo = e 2Arp °/ p r(0;2iVP /P)/ln2 



(14) 



at the laser threshold. This limit is "universal" , in the 
sense that it is independent of the degree of disorder 
(as long as we remain in the diffusive regime). We 
have Coo « P/2NP hx2 for P < NP Q and Coo « 
log 2 (P/2iVPo) ~ 7/ln2 for P > NP . The increase in 
the capacity by amplification in the diffusive regime is 
therefore up to a factor 3L/81 for P -C -ZVPq and up to 



a factor 1 + (biL/^QnP/NPo)- 1 for P > NP (L/l). 
All this is in contrast to the case of a waveguide with- 
out disorder, where the capacity is independent of the 
amplification. 

We now relax the requirement of heterodyne detec- 
tion and instead consider the maximum communication 
rate for any physically possible detection scheme || . We 
still assume that the information is encoded in coher- 
ent states, and use the same Gaussian a priori distri- 
bution p{n) oc exp(— \/i\ 2 Pq/P) as before. It has been 
conjectured |L8| that an input of coherent states with 
this Gaussian distribution actually maximizes the infor- 
mation rate for any method of non-entangled input with 
a fixed mean power (the so called one-shot unassisted 
classical capacity). 

The capacity for arbitrary detection scheme is given 
by the Holevo formula MM 



C H = H 



d 2 MP(»Pout(//) 



where H(p) = — Trplog 2 p is the von Neumann entropy. 
For a Gaussian density matrix p oc J d 2 /i exp(— |/i — 
pLo\ 2 fx) one has |21 



H(p) = (x + 1) log 2 (a; + 1) - x log 2 x = g(x), (15) 

Applying this formula to the Gaussian p out (/i) in Eq. (Q) 
and the Gaussian , we arrive at the capacity 



C n =g(rP/P + a-l)-g(a-l). 



(16) 



For a channel without amplification a — > 1 and so 
Ch = g( T P/Po), which lies above the capacity for hetero- 
dyne detection considered earlier. At the other extreme, 
upon approaching the laser threshold, a — > 00 and we 
have Ch — * log 2 (rP/(7Po), which is the same limiting 
expression as for heterodyne detection. 

The average over disorder can be carried out as be- 
fore by replacing a by a and averaging over t with the 
Rayleigh distribution V{f). The result is 



fP, a fP 
C H = — log 2 - — - + 

Po a — 1 Pq In 2 



1/a, 



-e^mil^ 1 )], (17) 



where TZ ef [/TZ' cS = 1 

As shown in Fig. ||, the dependence of Ch on the 
amount of amplification is non-monotonic — in contrast 
to the monotonically increasing C. Weak amplification 
reduces the capacity Ch, while stronger amplification 
causes Ch to rise again to the limit Coo at the laser 
threshold. The initial decrease for l a ^> L is described 
by 

Cn{L/l a ) pb C H (0) - (4lL 2 /3l 2 a ) log 2 (7Tl a /L). (18) 
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FIG. 3. Amplification dependence of the capacity C for 
heterodyne detection of coherent states [Eg. (pd|)] and the ca- 
pacity Ch for arbitrary detection [Eq. w/n)] . The input power 
is fixed at P/PqN = 1 and two values of l/L are chosen. 



investigating experimentally, for example in the context 
of a random laser . 

In conclusion, we have investigated the effect of ampli- 
fication on the information capacity of a disordered wave- 
guide, focusing on the competing effects of stimulated 
and spontaneous emission. We have compared the capac- 
ity C for heterodyne detection of coherent states with the 
Holevo bound Ch for arbitrary detection scheme. While 
amplification increases C for any magnitude of disorder 
and input power, the effect on Ch can be either favorable 
or not, as is illustrated by the "phase diagram" in Fig. |J. 

This research was supported by the Dutch Science 
Foundation NWO/FOM and by the RTN program of the 
European Community. 



Whether or not amplification ultimately increases Ch 
depends on the degree of disorder and on the input 
power. We indicate by A the region in parameter 
space where Coo > Ch(0) and by B the region where 
Coo < Ch(0). In region A strong amplification increases 
Ch while in region B it does not. The separatrix is 
plotted in Fig. 0. For P/NPq <C 1 the analytical ex- 
pression for this curve separating regions A and B is 
P/NPo = (3L/AI) exp(-3£/8Z+7), while for P/NP > I 
we find a saturation at l/L = 3/8e ps 0.14. This means 
that for P/NPq ^> 1 strong amplification increases the 
capacity Ch provided / < 0.14L. 
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FIG. 4. Curve in parameter space separating region A [in 
which Coo > Ch(0)] from region B [in which Coo < C*h(0)]. 
In region A amplification of sufficient strength increases the 
capacity Ch, while in region B it does not. 

At the laser threshold both C and Ch reach the same 
"universal" limit Coo given by Eq. (|l4|), which depends 
only on the dimensionless input power per mode P/NPq 
and not on the degree of disorder. This remarkably rich 
interplay of multiple scattering and amplification is worth 
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